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ABSTRACT 

Recent gamma-ray observations of the blazar lES 1101-232 (redshift z = 0.186) reveal 
that the unabsorbed TeV spectrum is hard, with spectral index a < 0.5 [F{y) oc 
We show that simple one-zone synchrotron self-Compton model can explain such 
hard spectra if we assume a power law energy distribution of the emitting electrons 
with a relatively high minimum energy. In this case the intrinsic TeV spectrum can 
be as hard as F{v) oc v^/'^^ while the predicted X-ray spectrum can still be much 
softer. The observations of lES 1101-232 can therefore be reconciled with relatively 
high intensities of the infrared background, even if some extreme background levels 
can indeed be excluded. We show that the other TeV sources (Mrk 421, Mrk 501 & 
PKS 2155-304) can be interpreted in the same framework, with a somewhat larger 
minimum energy. 

Key words: Radiative transfer - BL Lacertae objects: individual: Mrk 421, Mrk 501, 
PKS 2155-304, lES 1101-232 
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1 INTRODUCTION 

The intrinsic high-energy emission from extragalactic 
sources i s attenuated by the pair-production process 
jNishikovlll9"62 . Gould 1966, Steckc r et al. 1992 ). Photons 
in the TeV range (e-y) interact with low-energy photons 
from the extragalactic infrared background (em) and pro- 
duce electron-positron pairs (em + e^ — + e^e~). The to- 
tal absorption depends on the history of the reaction rate 
along the line-of-sight, and therefore, on the distance trav- 
elled and on the local density of the low-energy target pho- 
ton field, both a function of redshift z. This process will 
affect the high-energy TeV end of the observed blazar spec- 
tra. Since the optical depth depends strongly on the inci- 
dent TeV-photon energy the process leads to a steepening 
of the observed spectrum with respect to the intrinsic (un- 
absorbed) one. 

Thus, the observed spectra from distant blazars con- 
tain information on both the IR background radiation field 
history and the intrinsic properties of the source. From the 
viewpoint of blazar studies the latter is just an undesirable 
extrinsic effect that must be corrected for. In contrast, from 
the opposite perspective, it allows the measurement of the 
IR background. In both cases these two independent pro- 
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cesses have to be disentangled by either modeling both si- 
jriultaneously (e.g. iDe Ja gcr fc Steckcr 2002, Knciskc ct ^ 
j2p04), or measuring one independently fsee .Hauser fc Dwek 
|2p01 for a review). 

Unfortunately it proved very difficult to constrain the 
IR background independently, and different groups tried to 
constrain models for the IR background from TeV observa- 
tions assuming that the blazar sp ectrum is a well-behaved 
powe r-law (e.g. ICostamante et ali t2003. Dwek fc KrennrichI 
|2003)- The best perspective for such studies is offered by 
high-redshift blaza rs, sinc e the absorption will be strong. 
Recently, Aliaronia n et alj l|2006) presented a study along 
these lines, based on the TeV spectra of high-redshift blazars 
recently measured by the H.E.S.S. instrument. They de- 
reddened the observed TeV data of lES 1101-232 with dif- 
ferent models of the IR backgrounds, and concluded that 
the resulting intrinsic spectrum is always very hard. The in- 
dex of the intrinsic spectrum should be equal to a = 0.5 if 
the level of the IR background is close to the absolute lower 
limit obtained f rom direct integration of the light from re- 
solved galaxies jMadau fc Pozzettillioooll . Larger IR back- 
ground levels should correspond to intrinsic TeV spectra 
harder than a — 0.5. Note that this spectral index is already 
harder than the spectral index observed in other nearby 
TeV sources and that the diff erence will be mor e impo rtant 
for larger background levels. lAharonian et all i2006h also 
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commented on the theoretical difficulties to explain spectra 
harder than a — 0.5, which would imply quit e unexpected 
physics. These facts led lAharonian et al.l lj200(tl to prefer the 
solution of minimal IR background. 

In this paper we show instead that spectra harder than 
a = 0.5 can be obtained in "normal" one-zone homoge- 
neous synchrotron self-Compton (SSC) models, as long as 
there is a deficit of low energy electrons with respect to the 
extrapolation from higher energies. This can be achieved by 
assuming an emitting particle distribution which is a double 
power law, with a flat slope (< 2) of the first, low energy 
portion, or assuming a simple power law, but with a rela- 
tively high low energy cut-off. In these cases the limits on 
how hard a TeV spectrum can be is a = —1/3, while well 
accounting for a much softer X-ray spectrum. 

We discuss this result with the aim of finding limits on 
the IR background and also in the context of explaining the 
other TeV sources (with a softer intrinsic spectrum) with 
the same model. 
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2 SSC MODELS 

The simplest possible scenario able to explain the X-ray and 
TeV emission of a blazar assumes a homogeneous, spherical 
source filled by tangled magnetic field and relativistic elec- 
trons. 

Since in TeV blazars we must have electrons of very 
large energies, the scattering process occurs partly in the 
Klein-Nishina regime. In this regime the scattering process 
becomes inefficient, and for typical electron distributions the 
SSC process can be well described by approximating the 
scattering cross section with a step function, equal to the 
Thomson cross section for incident photon energies hu < 
rriaC? and zero otherwise^. This means that electrons of 
increasingly large energies scatter an increasingly smaller 
fraction of synchrotron seed photons, steepening the result- 
ing SSC spectrum with respect to the "canonical" value of 
Q = (n — l)/2. Note that this steepening introduced by the 
Klein-Nishina effects is virtually inescapable, since the scat- 
tering occurs completely in the Thomson limit only for ex- 
tremely large values of the beaming factor: & > 1000 (see e.g. 
[Katar zynski ct al. 2005a). In the next section we will find 
what is the high energy spectral index corresponding to a 
given particle energy distribution, assuming "usual" param- 
eters for our TeV source, namely: source radius R = IQ^^^^^ 
[cm], magnetic field intensity B = 1 — > 0.01 [G], Doppler 
factor 5 = 10 ^ 50 and approximating the particle en- 
ergy distribution by a broken power law with the index ni 
[A'^(7) = 7^7""^] below the break energy (7brk) and 712 above 
the break [N{"f) — K"f^^^"^'y~"^], where K is the density 
of the particles with 7 = 1. If the minimal 7 is greater than 
unity (7min > 1) the density of the particles with minimum 
energy [K') is given by K' = K-y~]^^ . 



The models shown in all figures and all the numerical results 
made use the full Klein-Nishina cross section; this approximation 
is mentioned here for illustrative purposes only. 
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Figure 1. Synchrotron emission (upper panel) and self-Compton 
radiation (lower panel) for different values of ni . The green spec- 
tra, generated by the truncated single power law particle en- 
ergy spectrum (ni = —00), provide the best fit for observa- 
tions of lES 1101-232 made by BeppoSAX (Wolter et al. 2002) 
and H.E.S.S. (Aharonian et al. 2006). Note that the spectra 
on the lower panel have been normalized to the peak level and 
the H.E.S.S. observations has been unabsorbed according to the 
model proposed by Kneiske et al. (2004). 



2.1 Flat electron distributions 

A fiat slope of the electron distribution (ni < 2) may help 
to explain the intrinsic TeV spectrum of l ES 1101-232. As- 
sumin g the lowest possible absorption level. I Aharonian et alJ 
(2006) have estimated the intrinsic spectral index a ~ 0.5 for 
this particular source. We verified that it is roughly possible 
to reproduce such a spectrum if we assume that ni ~ 1.5 
(note that in the Thomson limit a = 0.5 would correspond 
to ni = 2.). In this case, however, we do not obtain a power 
law intrinsic TeV spectrum. The spectrum is curved, with a 
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changing from a = 0.25 in the GeV range to a 
peak of the emission at a few TeV. 
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Since ni ~ 1.5 can explain the spectrum derived for 
the minimal IR background case, we then ask if a smaller 
value of m, producing an harder spectrum, can allow for 
a larger level of the IR background. In Fig. we show 
different synchrotron and inverse-Compton spectra for dif- 
ferent ni indices. The spectral index of the first part of 
the thin synchrotron emission (qi) depends directly on ni 
[ai = (ni — 1)/2] down to the limiting value m = 1/3. Below 
this limit the index remains constant ai = —1/3. This is due 
the fact that, below the limit, the dominant part of the low 
frequency emission is generated by electrons with 7 ~ 7brk. 
Therefore, the spectrum in the discussed range is equivalent 
to the spectrum produced by a monoenergetic population 
of electrons, which can be approximated by a power law 
with a = —1/3 below the maximum of the emission (e.g. 
iRvbicki &: LightmanlllQTsIl . This means that even in the ex- 
treme case of the truncated, single power law energy distri- 
bution (m = —00) there is always a "source" of low energy 
photons due to the low frequency emission (-Psyn ^ v ) 
of the high energy electrons (7 ~ 7brk) • The spectral index 
of the inverse-Compton emission in the Thompson regime 
is equivalent to the index of the target radiation field (e.g. 
lR,vbicki fc LightmanlllQTol) . Since ai is limited to -1/3 also 
the index of the inverse-Compton emission below the peak is 
limited to this value. Moreover, this index does not depend 
on 712. The index of the inverse-Compton emission above 
the peak is given by a = 2a2 — oi, where a\ — —1/3. This 
relationship can be easily derived from the appro ximation 
of the scattering derived bv lTavecchio et all l)l998h . 

To directly check the applicability of the SSC scenario 
that assumes ni < 2 we have applied our computations 
to the emission of lES 1101-232 observed by BepnoSAX 
l|Wolter et alJ l200(t and H.E.S.S. (|Aharonian et al.l \200(t ) 
experiments. The TeV data (here and in all the other cases 
in the following) have been unabsorbed according to the 
"best-fit" model proposed recently bv .Kneiskc et aL, 
The predictions of this mo del are also consisten t with those 
of the model proposed bv lStecker et alJ (120051) . The model 
provides a significantly stronger absorpti on level than the 
minim al possible absorption suggested bv lAharonian et alJ 
i2006h . However, our best fit for the expected intrinsic spec- 
trum, in a framework of the above mentioned SSC scenario, 
has been obtained for the limiting case (ni = —00). This 
indicates that the absorption is close to the highest pos- 
sible value that can be accommodated by this scenario. If 
the absorption is stronger , as suggested by some models (e.g. 
iDe Jaeer fc Stecke J200I1 . then our SSC scenario cannot ex- 
plain the intrinsic TeV emission of lES 1101-232. The de- 
tailed values of the physical parameters used in our compu- 
tations are discussed in the last section of this work, where 
we apply this scenario also to other sources. 

The truncated particle energy distribution corresponds 
to the limit ni — > —00 of a broken power law. For clarity, we 
will refer in the following to the break energy as 7min. For 
such a electron distribution, the shape of the TeV emission 
depends strongly on 7min, and we will therefore investigate 
this is some detail. 
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Figure 2. Synchrotron emission (upper panel) and the self— 
Compton radiation (lower panel) for different values of the 7min 
parameter. The best fit for the BeppoSAX (Pian et al. 1998) and 
CAT (Djannati-Atai et al. 1999) observations of Mrk 501 is indi- 
cated by the green lines. Note that the inverse-Compton spectra 
have been normalized and the TeV observations have been unab- 
sorbed according to the model proposed by Kneiske et al. (2004). 



2.2 Changing 7inin 

As we have shown in the previous subsection a particle dis- 
tribution with a low energy cut-off provides the best ex- 
planation for the emission of lES 1101-232. We here show 
the effects of changing the value of the low energy cut-off, 
which we will call 7miii. The result of our calculations is 
s hown in Fig.|21 We have applied this test for the BeppoSAX 
(Pian et aljil99al and CAT (Diannati-A tai et al. 1999) ob- 
servations of Mrk 501 from 16"^ April, 1997. 

For this test, we fix the slope of the particle distribu- 
tion to n = 2, which provides the best fit for the used X-ray 
observations of Mrk 501. In the first approach we assume 
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Figure 3. X-ray and TeV emission of four TeV blazars located at different redshifts. For all the sources we have applied the same SSC 
scenario, assuming that the emitting electrons are distributed in energy as a single power law with a low energy cut-off at 7inin- AH 
the TeV spectra have been unabsorbed according to the model of Kneiske et al. (2004). Solid lines are the best fits for the expected 
intrinsic emission whereas dashed lines indicate the observed spectra. For an easier comparison of the different sources, we transformed 
the observed flux in luminosity. 



the minimal possible value of 7min = 1- This means that 
the model is equivalent to the standard SSC scenario and is 
providing a relatively broad peak for the inverse-Compton 
emission. Nevertheless, this approach is able to roughly ex- 
plain the unabsorbed TeV spectrum. However, the increase 
of the 7miii value by about three orders of magnitude pro- 
vides a slightly narrower peak and therefore better explana- 
tion for the spectrum. A value of 7min greater by one more 
order of magnitude is the limiting value if we want to explain 



the entire X-ray spectrum. However, this already high value 
of 7min is drastically modifying the TeV spectrum, where 
the limiting slope (q = —1/3) of the emission becomes well 
visible (blue line in Fig. 2). Finally, a quasi-monoenergetic 
electron distribution cannot explain the X-ray nor the TeV 
emission. Note that for values of 7max close to 7min, the sec- 
ond part of the TeV spectrum, above the peak, is curved 
and it is not a power law. 
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3 APPLICATION TO OTHER TEV BLAZARS 

We have already applied our specific SSC scenario to 
lES 1101-232 and Mrk 501. In this section we apply this 
model to two other TeV sources, Mrk 421 and PKS 2155- 
304, and compare the results. 

The RXTE and HE GRA observations of Mrk 421 
iKrawczvnski et al.ll200]l) can be fitted assuming that the 
X-ray emission is done by the exponential synchrotron tail 
of the highest energy particles at 7max. Therefore, the index 
of the particle energy spectrum and 7min are not constrained 
by the data. The same approach can explai n the RXTE and 
H.E.S .S. observations of PKS 2155-304 jAharonian et alJ 
[200^. Note that it is also possible to explain these obser- 
vations, as well as the observations of Mrk 501, assuming 
7inin — 1. We performed our computations to show that 
the assumption of having 7min 3> 1 (required to explain 
the emission for lES 1101-232), works very well also for 
other sources. In Fig.|21we show the results of our modeling 
and give the values of the input parameters. Note that for 
Mrk 501 we provide also the spectral fits for the observations 
made on April 7^^, 1997. Moreover, in the case of Mrk 421, 
Mrk 501 and PKS 2155-304 the X-ray and TeV observations 
were made almost simultaneously, whereas for lES 1101-232 
we use the X-ray spectrum from January 1997 and the TeV 
data collected over March- June 2004-2005. 

Note that the value of the magnetic field used in our 
calculations is significantly smaller that the value estimated 
from the equipartition between the magnetic field energy 
density and the electron energy density. 



4 DISCUSSION 

We have shown that the standard SSC model can explain 
the specific observations of the relatively distant TeV blazar 
lES 1101-232. Assuming a truncated power law energy dis- 
tribution, with 7min ^ 1, we can well explain the intrinsic 
TeV spectrum even assum ing a stronger absorptio n than the 
upper limit suggested bv lAharonian et al.l ll2006f) . This as- 
sumption appears in very good agreement with the recent 
estimates of the IR interg alactic background iKneiske et alJ 
|2004 IStecker et aklliooai . The limit on how hard an SSC 
spectrum can be is a — —1/3. This limit in turn translates 
in a limit on the amount of absorption sufi'ered by TeV pho- 
tons, and hence on the level of the IR background. 

Concerning the reason of having 7min 2> 1, we can pro- 
pose two answers. The first is the injection of an already 
truncated power law spectrum from a shock region. If, in ad- 
dition, the cooling process is inefficient for the the particles 
at lower energies (close to 7min), the particle distribution 
will not develop a low energy tail in a "reasonable" time, 
namely a dynamical time, over which we can assume that 
the physical conditions (size, densities, magnetic field) are 
constant (since the source has bulk motion, after a dynam- 
ical time it will be located in a larger part of the jet, and 
not contribute any longer to the observed flux). Note that 
in our particular simulations also particles with 7 ~ 7max 
are not affected by the cooling in one crossing time (R/c). 

In a second, more complex, scenario all particles are ac- 
celerated efficiently at the shock front almost up to some 
maximal energy. Then they escape downstream of the shock 



and cool. However, the cooling process may be compen- 
sated (at l east in part) by some st ochastic turbulent ac- 
celeration jKatarzvnski et al]l2005bh . The competition be- 
tween re-acceleration (heating) and cooling may single out 
an equilibrium energy where the two processes balance. This 
equilibrium energy may correspond to the 7min assumed in 
our simple SSC scenario. 

The changes of 7min and 7max may have various impacts 
for the source variability, especially for the correlation be- 
tween the X-ray and TeV emission. As an example, when 
722 < 3 and 7min ^ 7max the change of 7inin may not affect 
the X-ray emission while having a strong impact for the 
TeV radiation. Such a scen ario could explain the or phan 
TeV flare in lES 1959-f650 iKrawczvnski et ahlliooi) . On 
the contrary, a change of 7max would produce a significant 
X-ray variability without any change in the TeV flux below 
the TeV peak. 



ACKNOWLEDGMENTS 

We are grateful to A. Djannati-Atai, H. Krawczynski and 
B. Giebels for the data obtained by the CAT, HEGRA and 
H.E.S.S. experiments. We acknowledge the EC funding un- 
der contract HPRCN-CT-2002-00321 (ENIGMA network). 
JG thanks for the hospitality at the Osservatorio di Brera. 



REFERENCES 

Aharonian F., Akhperjanian A.G., Bazer-Bachi A.R., et 

al. 2005, A&A, 442, 895 
Aharonian F., Akhperjanian A.G., Bazer-Bachi A.R., et 

al. 2006, Nat, in press (astro-ph/0508073 I 
Costamante L., Aharonian F., Ghisellini G. & Horns D. 

2003, New Astron. Rev., 47, 677 

De Jager O.C., & Stecker F.W. 2002, 566, 738 
Djannati-Atai A., Piron F., Barrau A., et al., 1999, A&A, 
350, 17 

Dwek E. & Krennrich F., 2005, ApJ, 618, 657 
Gould R. J., Schreder G. P., 1966, Phys. Rev. Lett. 16, 252 
Hauser M. G., Dwek E., 2001, ARA&A, 39, 249 
Katarzynski K., Ghisellini G., Tavecchio F., Maraschi L., 

Fossati G. & Mastichiadis A. 2005a, A&A, 433, 479 
Katarzynski K., Ghisellini G., Mastichiadis A., Tavecchio 

F. & Maraschi L. 2005b, A&A, in press 
Kneiske T.M., Bretz T., Mannheim K., & Hartmann D. H., 

2004, A&A, 413, 807 

Krawczynski H., Sambruna R., Kohnle A., et al., 2001, 
ApJ, 559, 187 

Krawczynski H., Hughes S.B., Horan D., et al. 2004, ApJ, 
601, 151 

Madau P. & Pozzetti L., 2000, MNRAS, 312, L9 
Nishikov A. I., 1962, Sov. Phys. JETP, 14, 393 
Plan E., Vacanti G., Taghaferri G., et al., 1998, ApJ , 492, 
L17 

Rybicki G. & Lightman A. P., 1979, Radiative Processes in 
Astrophysics, Wiley Interscience, New York 

Stecker F. W., de Jager O.C. & Salamon M. H., 1992, ApJ, 
390, L49 

Stecker F. W Malkan, M. A. & Scully, S. T., 2005, ApJ, 
submitted, ( |astro-ph/05 1 0449 } 



6 K. Katarzynski et al. 

Tavecchio P., Maraschi L. & GhisoUirii G., 1998, ApJ , 509, 
608 

Wolter A., Tavecchio F., Caccianiga A. & Tagliaferri G., 
2000, A&A, 357, 429 



